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We compute the entropic interactions between two colloidal spheres immersed in a dilute suspen- 
sion of semi-flexible rods. Our model treats the semi-flexible rod as a bent rod at fixed angle, set by 
the rod contour and persistence lengths. The entropic forces arising from this additional rotational 
degree of freedom are captured quantitatively by the model, and account for observations at short 
range in a recent experiment. Global fits to the interaction potential data suggest the persistence 
length of fd- virus is about two to three times smaller than the commonly used value of 2.2 jj,m. 
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Colloidal dispersions exhibit a fascinating range of 
equilibrium and non-equilibrium structures, and they 
have important impact on our daily lives The inter- 
actions between suspension constituents determines the 
stability of the dispersion against fiocculation, and the 
phase behavior of the colloid. Quantitative models and 
measurements of these interactions test our basic under- 
standing about these systems, and enable experimenters 
to better control suspension behaviors and properties. In 
this paper we focus on a particular class of entropic inter- 
action, exploring the forces between spheres in a suspen- 
sion of rodlike particles. This system class has produced 
a variety of interesting phases [2, pL m, and has stimulated 
several theoretical models || 6, 0, [| |t| and a measure- 
ment jl^ of the rod-induced depletion interaction. 

The depletion attraction between two spheres im- 
mersed in a dilute suspension of thin rods of length, L c , 
was first considered by Asakara and Oosawa juT. Their 
most important physical insight was that rods in sus- 
pension gain both translational and rotational entropy 
when the sphere surfaces come within L c of one another. 
Subsequent theories computed the attraction more ac- 
curately within the Derjaguin approximation j^, [?j and 
beyond [pi]. However, in many practical scenarios the 
rods are not rigid, and current theories do not account 
for the flexibility of the rods. Indeed, flexibility effects 
can be important as evidenced by a recent interaction 
measurement 10 of micron diameter spheres in suspen- 
sions of fd-virus; in this case systematic deviations be- 
tween experiment and "rigid-rod" theories were found at 
short-range, and were suggested to arise as a result of 
the flexibility of the fd-virus. Flexible or bent rods have 
an additional degree of freedom: the rotation about their 
central axis. As the spheres get closer, this degree of free- 
dom is depleted, the system entropy increases, and the 
sphere interactions become even more attractive. 

A quantitative model for this observation is still lack- 
ing, and indeed a complete theory of semi-flexible rods 
near surfaces remains a difficult task. In this paper, we 
introduce a simple model to compute the depletion po- 
tential between two spheres in a dilute solution of semi- 
flexible rods. We use the model to quantitatively explain 
the experiments of Ref. p0[. The model accounts for 



the entropic effects of flexibility at short-range, and pro- 
vides an accurate fit of the measured interaction poten- 
tials. The model also provides a means to extract the 
persistence length l v and the contour length L c of the 
suspended semi-flexible rods from interaction potential 
data. Global fitting of the data suggests that the per- 
sistence length of fd-virus is two to three times smaller 
than the commonly used value of 2.2 /im fll2| ]. 

Our model relies on the assumption that if the rods are 
sufficiently stiff, they may be accurately approximated by 
two rods of length L = L c /2, attached together at a fixed 
angle tt — 2a, as shown in Fig. m. The angle a may be 
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estimated by a = cos -1 R/L c , where R = ((R 2 )) is 
the average end-to-end distance. R is related to l v and 
L c by P 
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This approach simplifies the problem, while still captur- 
ing the essential physics. In particular, we show that the 
part of the depletion potential associated with new rota- 
tional degrees of freedom is short ranged, i.e. of order 




FIG. 1: (a) Typical configurations of a semi-flexible rod whose 
contour length L c is comparable to the persistence length l v . 
(b) Our approximation of semi- flexible rod in (a); two stiff 
rods of length L attached at a fix angle, (c) A bent rod near 
a flat wall and (d) confined between two walls. 
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the bent rod width, W = Lsina. Importantly, when W 
is significantly less than the particle diameter, the rota- 
tional part of the depletion interaction can be treated 
within the Derjaguin approximation J]J]. 

In the presence of repulsive walls (see Fig. [I]), the 
rotational degrees of freedom of a bent rod are restricted. 
Consider a bent rod with one end displaced by z from 
the wall and with orientation (u, <fr). The probability 
of finding a rod in such a configuration is given by the 
Boltzmann factor: /(r, u, <j)) oc exp[— (3U e (r, u, </>)], where 
[3 = 1/kgT, kg is the Boltzmann constant, and T is the 
temperature. The hard wall potential U e is infinite if 
any part of the rod touches the wall and is otherwise 
zero. We consider the case where the concentration of 
the rods is sufficiently low so that the thermodynamics 
are well characterized by the Grand potential of an ideal 
gas of rods 
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Here N is the number of rods, 
tension by the difference 
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Here po is the average density of the rods and S is the 
surface area of the wall. To compute the integral in Eq. 
(||) , we enumerate all the configurations of the bent rod 
just touching the walls. 

Let us first consider a single flat wall, as shown in Fig. 
[l]. There are three regions to consider: (i) < a < tt/A, 
(ii) 7r/4 < a < 7r/3, and (iii) ir/3 < a < tt/2. For 
a < 7r/4, we observe that when 8\{z,a) < 9 < 9 2 {z 1 a), 
for 



9i(z,a) = cos 1 
%{z, a) = a + cos~ i z/L, 



2L cos a 
-1 



(4) 
(5) 



the rotation of the rod about its symmetry axis is re- 
stricted to A a < 6 < 2tt — d> a , where 



b a (z,9) = cos 1 



L sin a sin 9 



cot 9 cot a 
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FIG. 2: The surface tension of a bent rod in the presence of 
a flat plane wall. 



Using this construction, the surface tension is 
f3Aj(a) Lcosa 
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dx 



d9 sin 9 (j) a (x, 9), 
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where x — z/L. Similarly, for ^ > a > j, we have 
f3Aj(a) Lcosa L 
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a — cos 1 x 



In Fig. H we plot A7(a) as a function of a. Note that 
the limiting values, A7(0) = ^poksTL and A7(-7r/2) = 
jPoksTL, agree with previous results fll4| . 

We now turn to the calculation for two walls (see Fig. 
1J). Since the rods are stiff and a is small in the exper- 
iment of interest, we focus on the case where a < 7r/4. 
For a given separation of the walls d, we divide the inter- 
val < z < d/2 into different regions, wherein 61,82,03, 
and #4 take on different values. The new angles are 
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FIG. 3: The depletion interaction V(d) [Eq. ([TJ)] between 
two planar walls (the solid curve) mediated by a bent rod of 
contour length L c — 2L with (a) a — tt/7 and (b) n/5. The 
dashed curve is the depletion interaction of a straight rod with 
R = 2Lcosa. At large distances, they show little difference 
but the restriction on the additional degree of freedom at 
shorter distances gives rise to a stronger attraction in V(d), 
which is bounded below by the potential of a straight rod 
of R = 2L (the dotted line). This is qualitatively the effect 
observed in experiment of Ref. (See also Fig. §). 



3 



9i{z, a) = n — cos 



2L cos a 



(10) 



If 6x < 9 < 02 and 82 < 6*3, 4> is restricted to 4> a < <j> < 
2tt — <f) a . If 9 3 < 9 < 84, (j) is restricted to < (j) < ^ — 4>b 
and tt + (f>f, < 4> < 2tt with 



<j>b(z, 9) = cos 1 



1 



L cos(7r — 9 — a) — (d — z) 



Lsinasin(7r — 9) 

(11) 

When 92 > #3, 4> is further restricted to <j) a < 4> < tt — 
& and 7r + 06 < < 27r — a if 6» 3 < 6> < 6* 2 . Thus, 
the depletion potential per unit area defined by V(d) = 
A7[d, a] — A7[oo, a] is 
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Here ' indicates integrations over phase space restricted 
to the allowed values. Fig. 3 depicts the depletion poten- 
tial between two walls for different a. At large distances 
the potential is determined by the "end-to-end" distance 
R of the rod. At short distances the rotational degree of 
freedom becomes important and increases the attraction 
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FIG. 4: Interaction potential between a pair of 1.0 /Km sil- 
ica spheres in a suspension of fd virus with concentration 
0.67mg/ml. The dotted (dashed) lines are generated by the 
YJM model with R = L c = 900 nm (R = 740 nm; L c = 900 
nm) . The solid lines are generated by Eq. ( fl4| ) with R = 740 
nm and L c — 900 nm. Clearly, the agreement of experimen- 
tal data and our model which includes the addition rotational 
degree of freedom of a bent rod is excellent. The dash-dotted 
vertical line indicates W = 0.23 fim. 



between walls. The potential is bounded below by the 
potential of a straight rod with length 2L. Although our 
calculation has been done for two walls, we expect the 
same qualitative features to hold for two spheres. 

V (d) can be written as a sum of 2 pieces. The first term 
is the depiction potential of a straight rod with length 
R = 2L cos a jL4| . The second term depends only on the 
additional rotational degree of freedom of the bent rod. 
Moreover, the range of T(d, a) is of order of the width of 
the bent rod, W = Lsina, which is small compared to 
the sphere radius in Ref. [jl(|. These observations sug- 
gest that to approximate the depletion potential for two 
spheres, the latter term may be treated in the Derjaguin 
approximation, while the first term replaced by the YJM 
rigid-rod model |B. Thus, we write 
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where a is the sphere radius and h their closest surface 
separation. K(h/R;a/R) is the potential between two 
spheres due to a straight rod of length R, which reduces 
to the Derjaguin expression Krj(h/R) = -|(1 — h/R) 3 in 
the limit a/R » 1 |J. 

Fig. |i| displays a typical experimental data set of Ref. 
|]i"of with three different models - (i) the YJM model 
(dotted line), whose potential is given by the first term 
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FIG. 5: Interaction potential between pairs of (a) 1.0 /im and 
(b) 1.6 (J,m silica spheres in a suspension of fd virus with differ- 
ent concentration. The solid lines are generated by Eq. (|l4|), 
with best fit parameters that give smallest x 2 - 
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in Eq. (14) with R = 900 nm, the contour length of fd, 
(ii) the YJM-KP model (dashed line), whose potential is 
given by the first term in Eq. (14) with R = 740 nm, 
and (iii) the bent rod model (solid line), U s (h) in Eq. 
@ with R = 740 nm and L c = 900 nm. The circles are 
experimental data for 1.0 /mi diameter silica particle in 
a dilute (0.67mg/ml) solution of fd virus. The theory 
curves are computed with no free parameters and are 
then numerically blurred to account for the instrument's 
spatial resolution (see Ref. [f[o)). Clearly, our model gives 
the best fit to the experimental data. In particular, both 
YJM and YJM-KP models, while having approximately 
the right magnitude and shape, fail to account for the 
overall curvatures of the experimental curve. Further, 
while the YJM-KP model agrees with most of the data 
at large h, our model clearly accounts for the depth of 
the measured potential near contact. 

In order to explore the best fits more quantitatively, 
we computed the x 2 value of our models for all data sets 
with R ranging from 720 - 825 nm and L c = 880, 900, 
and 920 nm. If a fixed concentration (measured experi- 
mentally) is assumed, x 2 is smallest for R = 780 nm and 
L c = 920 nm. If the concentration is allowed to vary 
within its ±5% experimental error, then % 2 is small- 
est for R = 740 nm and L c = 900 nm. Fig. |5| shows 
best fits for each of concentration. Note that the width 
W ~ 200 nm, is smaller than the radius of the colloidal 
spheres. This justifies a posteriori the Derjaguin approx- 
imation made in Eq. (|l4|). Furthermore, we can estimate 
t p using Eq. (|l|) and the values for R and L c above, yield- 



ing £ p ~ 850 nm (fixed concentration) and £ p ~ 680 nm 
(variable concentration). Our results for L c of fd are con- 
sistent with the literature, i.e. 850 nm < L c < 920 nm 
|]l5| . However, our values for £ p should be contrasted to 
the often-quoted value, £ p = 2.2 /im Jl|, |l7). The lat- 
ter is based on a fitting of dynamic light scattering data 
with theoretical models |l7], [ll| , whose assumptions may 
well be questioned in the light of our results. Indeed, 
smaller values of £ p have also been reported based on dy- 
namic structure factor models of semi-flexible filaments 
]l9[ , and using electron microscopy pofl . 

We have presented a simple analytical model for the 
depletion interaction between two spheres mediated by 
semi-flexible rods, and demonstrated its quantitative 
agreement with experimental data. Our theoretical 
model combined with interaction measurements provides 
a basis for extracting the persistence length of a semi- 
flexible rod. 
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